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ABSTRACT 


A  preliminary  design  study  was  condueted  to  determine  the  geometry  for  an 
enhaneed  mixing  eduetor  system  for  the  Landing,  Helicopter  Assault  Ship  Replacement 
(LHA  (R))  program  gas  turbine  exhaust.  A  one-dimensional  analytical  model,  with  a 
correction  factor  applied  to  the  secondary  mass  flow,  was  developed  to  predict  the 
secondary  air  mass  flow  rate  and  the  exhaust  temperature  at  the  mixing  tube  exit  plane. 
The  resultant  design  consisted  of  a  high  aspect  ratio  lobed  nozzle  and  a  mixing  tube.  The 
model  was  also  used  to  predict  the  backpressure  developed  by  the  ducting  configuration. 
The  proposed  design  resulted  in  a  50%  reduction  in  exhaust  temperature  with  only  a  6 
inch  H2O  increase  in  backpressure.  A  detailed  design  of  the  oval-to-rectangular 
transition  duct  is  provided,  based  on  empirical  data  from  a  similar  duct  design.  The  study 
also  included  a  prediction  of  plume  radiation  intensity  in  the  3-5  pm  band  for  various 
aspect  ratio  slots. 
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I.  INTRODUCTION 


A,  MOTIVATION 

The  search  for  a  more  efficient  propulsion  plant  that  meets  or  exceeds  current 
capability,  while  reducing  maintenance  requirements  has  led  the  United  States  Navy  to 
investigate  an  all-electric  ship.  A  first  step  on  the  road  to  the  all-electric  ship  is  a 
modification  to  the  WASP  class  Landing,  Helicopter  and  Docking  Assault  ship  (LHD), 
Reference  [1].  The  eighth  ship  of  the  WASP  class  (LHD  8)  is  to  incorporate  a  hybrid  gas 
turbine  and  electric  motor  for  propulsion  into  the  existing  hull  and  superstructure  of  the 
LHD  1 .  The  WASP  class  (LHD  1  through  LHD  7)  propulsion  plant  consists  of  two  600 
psi  boilers  and  two  35,000  hp  steam  turbines  that  can  propel  the  40,500  ton  Amphibious 
Assault  Ship  at  25  knots.  In  order  to  meet  the  Horsepower  requirements  of  the  LHD  8, 
the  General  Electric  LM2500+  was  chosen  to  replace  the  steam  plant  along  with  a  5,000 
hp  auxiliary  propulsion  motor  that  will  be  used  for  loitering  on  station.  The  Landing, 
Helicopter  Assault  Ship  Replacement  (LHA  (R))  program  will  be  the  follow-on  to  the 
aging  WASP  class,  and  is  to  have  the  same  propulsion  plant  and  similar  hull  form  as  the 
LHD  8  with  improvements  in  superstructure  design  in  order  to  reduce  radar  cross  section. 

There  are  many  advantages  to  using  gas  turbines  for  propulsion.  Some 
advantages  include,  improved  power-to-weight  ratio,  reduced  maintenance  requirements, 
and  reduced  engineering  watch  stander  requirements.  However,  the  temperature  of  the 
gas  turbine  exhaust  can  be  in  excess  of  1000  °F  (800  K),  roughly  500  °F  (530  K)  higher 
than  a  conventional  steam  plant.  The  increased  temperature  can  damage  electronic 
equipment  mounted  on  the  mast  as  the  hot  exhaust  plume  comes  in  contact  with  it. 
Additionally,  the  ship  surfaces  which  normally  radiate  as  nearly  a  black  body  in  the  8-14 
pm  band  are  heated  after  coming  in  contact  with  the  plume,  and  now  radiate  at  a  higher 
intensity  given  by  I=(saT"^)/7i  (where  T  is  the  exhaust  temperature  in  Kelvin,  s  is  the 

o 

emissivity  and  a=5. 67x10'  is  the  Stefan-Boltzmann  constant).  Figure  [1]. 

Since  the  intensity  of  the  plume  also  goes  as  T"^,  the  intensity  of  gas  turbine 

exhaust  can  be  16  times  greater  than  the  conventional  steam  plant.  The  two  main 

byproducts  of  combustion,  CO2  and  H2O,  strongly  emit  spectral  radiation  in  the  3  -  5  pm 

band,  and  though  the  atmosphere  is  a  good  absorber  of  radiation,  it  has  gaps  in  the  3-5 
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|im  and  8-12  |a,m  bands.  Consequently,  the  addition  of  a  gas  turbine  power  plant  to  U.  S. 
Navy  warships  inereases  the  ships  suseeptibility  to  deteetion,  tracking,  and  targeting  by 
infrared  seeking  anti-ship  missiles.  It  is  possible  to  redesign  the  ship  so  that  the  mast  is 
far  enough  away  from  the  hot  exhaust  plume  that  damage  to  electronics  does  not  occur. 
However,  simply  moving  the  mast  in  relation  to  plume  does  not  address  the  increased 
infrared  signature  of  the  ship.  Developing  a  means  to  reduce  the  gas  turbine  exhaust 
temperature  will  address  both  problem  of  damage  to  equipment  and  increase  infrared 
signature. 

Northrop  Grumman  Ship  Systems  (NGSS)  has  provided  a  design  of  an  active 
cooling  system  for  LHD  8  that  consists  of  two  150  hp  axial  fans  that  are  each  60  inches  in 
diameter.  The  system  proposed  by  NGSS  requires  a  significant  amount  of  space,  which 
results  in  increased  weight,  complexity  and  maintenance  requirements.  The  addition  of 
the  fan  room  near  and  underway  replenishment  (UNREP)  station  reduced  the  amount  of 
space  available  for  cargo  handling,  and  also  generated  an  unacceptable  noise  hazard.  To 
avoid  the  noise  hazard,  reduction  in  space  near  the  UNREP  station,  and  significant 
increase  in  topside  weight,  it  was  decided  to  design  a  passive  cooling  system. 
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Figure  1  Black  Body  Emissive  Power,  from  Reference  [14], 

B,  BACKGROUND 

The  challenges  posed  in  cooling  gas  turbine  exhaust  are  not  new.  Wilsted, 
Huddleston  and  Ellis  [2]  of  the  National  Advisory  Committee  for  Aeronautics  (NACA), 
now  known  as  the  National  Aeronautics  and  Space  Administration  (NASA),  researched 
the  effect  of  temperature  on  the  performance  of  ejectors  for  gas  turbine  exhaust  in  1949. 
In  1953,  Hussmann  [3]  prepared  a  report  for  the  United  States  Navy  Bureau  of  Ships, 
now  known  as  the  Naval  Sea  Systems  Command  (NAVSEA),  that  detailed  how  to  design 
an  eductor  system  for  gas  turbine  powered  ships.  In  1975,  the  Navy  commissioned  USS 
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SPRUANCE  (DD  963),  the  first  large  U.  S.  warship  to  use  gas  turbines  for  propulsion. 
The  gas  turbine  exhaust  eooling  system  employed  on  SPRUANCE  elass  destroyers  is  a 
single  nozzle  eduetor,  similar  to  the  sehematic  shown  in  Eigure  [2],  and  a  Boundary 
Eayer  Infrared  Signature  Suppression  (BEISS)  eap.  The  eduetor  mixes  the  hot  exhaust 
gases  (1)  with  eool  ambient  air  (2),  also  known  as  primary  flow  and  seeondary  flow 
respeetively.  The  BEISS  eap  introduees  a  thin  layer  of  eool  air  between  the  exiting  flow 
and  staek  surfaee.  The  eooler  boundary  layer  reduees  the  surfaee  temperature  of  the 
staek  exit  and  provides  some  additional  eooling  of  the  exhaust  gases.  The  design  of  the 
eduetor  system  used  on  Navy  ships  has  improved  over  several  ship  programs.  Researeh 
at  the  Naval  Postgraduate  Sehool,  in  eonjunetion  with  David  Taylor  Researeh  Center 
(Annapolis,  MD)  and  NAVSEA  began  in  1977,  Ellin  [4],  to  enhanee  the  mixing  of  the 
hot  exhaust  gases  with  seeondary  eooling  air  by  using  multiple  nozzles  viee  a  single 
nozzle  in  the  eduetor  system.  The  result  of  this  researeh  is  the  eurrent  gas  turbine  exhaust 
eooling  system  used  on  the  TICONDEROGA  (CG  47)  Class  Cruisers  and  ARLEIGH 
BURKE  (DDG  51)  Class  Destroyers,  Eigure  [3]. 
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(3) 


Figure  2  Schematic  of  a  simple  eductor,  from  Reference  [6] 
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Figure  3  (a)  DDG  5 1  eductor  system,  (b)  Plan  view  of  DDG  5 1  eductor  nozzle,  after 
Reference  [9]. 
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Although  the  eductor  systems  currently  in  use  on  U.  S.  Navy  ships  are  able  to 
reduce  the  gas  turbine  exhaust  temperature  to  a  level  comparable  to  that  of  its  steam 
counter  part,  the  length  of  the  mixing  tube  is  greater  than  twice  the  diameter  of  the 
mixing  tube.  In  the  case  of  the  DDG  51,  this  results  in  a  mixing  tube  length  of  roughly 
20  feet.  Another  successful  eductor  system  developed  in  Canada  is  the  DRES  Ball, 
Figure  [4].  The  design  allows  ambient  air  to  be  drawn  into  the  core  of  the  hot  primary 
air,  and  slightly  reduces  the  amount  of  topside  weight. 


EXHAUST 

FLOW 


Figure  4  DRFS  Ball  device,  from  Reference  [15]. 

In  an  effort  to  reduce  the  mixing  tube  length,  and  thus,  topside  structure  and 
weight,  a  significant  amount  of  research  has  been  conducted  in  the  field  of  enhanced 
mixing  using  lobed  nozzles.  Figure  [5].  In  1987,  Presz,  Blinn  and  Morin  [5] 
demonstrated  that  use  of  a  lobed  nozzle  could  reduce  the  mixing  tube  length  to  diameter 
ratio  (F/D)  to  nearly  unity.  Fobed  nozzles  are  used  extensively  in  aircraft  applications 
due  to  their  short  mixing  length. 

Since  eductors  have  proven  to  reduce  the  gas  turbine  exhaust  temperature  to 
acceptable  levels,  it  was  decided  to  employ  the  same  technology  to  the  FHA  (R). 
Additionally,  a  lobed  mixer  will  be  used  for  the  nozzle  to  minimize  the  amount  of  topside 
space  and  weight  added  to  the  ship  from  using  an  eductor  system. 
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Figure  5  Schematic  of  a  lobed  mixer  illustrating  streamwise  vortices,  from  Reference  [5]. 

C.  OBJECTIVES 

The  objectives  of  the  research  are  to; 

1 .  Provide  a  preliminary  design  of  a  passive  cooling  system  for  the  LHA  (R)  gas 
turbine  exhaust  that  meets  or  exceeds  the  reduction  in  temperature  that  is 
provided  by  the  DDG  5 1  exhaust  cooling  system,  and  fits  within  the  existing 
geometry  of  the  LHD  8  uptake  space. 

2.  Estimate  the  total  pressure  loss  in  the  proposed  system  design  and  compare 
with  the  total  pressure  loss  in  an  uncooled  system.  The  proposed  design  will 
include  an  estimate  of  the  total  pressure  loss  a  circular-to-rectangular 
transition  duct. 

3 .  Study  the  effect  of  exhaust  plume  aspect  ratio  on  the  radiant  intensity  for  the 
DDG  51  gas  turbine  exhaust  system.  The  study  will  determine  if  there  is  a 
possibility  of  designing  the  system  with  a  significant  reduction  in  radiant 
intensity. 
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D,  ORGANIZATION 

Chapter  II  contains  the  I-D  model  development  for  the  eductor  system.  The 
discussion  will  include  assumptions  that  were  made  in  order  to  simplify  the  calculations. 
This  chapter  also  contains  a  brief  discussion  of  how  the  eductor  model  was  validated. 
The  pumping  ratio  correction  factor  is  explained  in  this  chapter. 

Chapter  III  contains  the  design  studies.  The  first  study  varied  the  mixing  tube 
cross-sectional  area,  and  observed  its  effect  on  backpressure.  The  second  study  varied 
the  nozzle  cross-sectional  area,  with  constant  mixing  tube  exit  area,  and  observed  its 
effect  on  exit  temperature  and  backpressure.  The  results  from  the  two  studies  were  used 
to  determine  the  design  point  for  the  eductor  system.  A  comparison  of  the  backpressures 
and  mixing  tube  exit  temperatures  for  the  proposed  design  and  an  uncooled  system  is 
contained  in  this  chapter. 

Chapter  IV  contains  the  transition  duct  design.  The  discussion  includes 
development  of  the  geometry  and  an  estimate  of  the  total  pressure  loss  associated  with 
the  duct.  Validation  of  the  code  used  to  determine  the  losses  is  included. 

Chapter  V  contains  a  study  of  how  changing  the  circular  exit  cross-section  of  the 
DDG  5 1  exhaust  system  to  a  rectangular  cross-section,  and  then  varying  the  aspect  ratio 
affects  the  radiant  intensity  of  the  exhaust  gases. 

Chapter  VI  contains  the  summary  and  conclusions. 
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II.  MODEL  DEVELOPMENT  AND  VALIDATION 


A.  OVERVIEW  AND  ASSUMPTIONS 

A  one-dimensional  model  of  an  ideal  eduetor  system  was  developed  to  determine 
the  basie  geometry  required  to  meet  the  objeetives  of  this  study.  Simulating  a  real,  three- 
dimensional  eduetor  system  with  a  1-D  ideal  model  required  that  several  assumptions  be 
made  ineluding,  the  flow  is  uniform,  one-dimensional,  and  uni-axial.  Additionally, 
frietion  on  the  mixing  tube  wall  was  negleeted.  Negleeting  mixing  tube  wall  frietion  and 
assuming  uniform  flow  will  overestimate  the  secondary  mass  flow  rate,  and 
underestimate  the  mixing  tube  exit  temperature.  To  account  for  the  errors  induced  in  the 
ideal  model,  a  correction  factor  that  reduces  the  predicted  entrainment  of  secondary  mass 
flow  will  be  introduced. 

The  following  assumptions  were  made  in  the  development  of  the  1-D  ideal 
eductor  model: 

1.  The  flow  is  uniform,  uni-axial  and  one-dimensional  at  the  nozzle  exit,  the 
mixing  tube  entrance  and  the  mixing  tube  exit. 

2.  The  flow  is  incompressible  since  the  maximum  velocity  in  the  system,  at  the 
nozzle  exit,  is  less  than  Mach  0.3. 

3.  The  exhaust  gases  are  ideal  gases. 

4.  The  properties  (density,  viscosity,  specific  heat,  etc...)  of  the  exhaust  gases 
are  that  of  air. 

5.  The  density  of  the  exhaust  gases  in  the  mixing  tube  varied  axially  due  to  the 
temperature  change  from  entrance  to  exit. 

6.  The  walls  of  the  exhaust  ducting,  including  the  mixing  tube,  are  adiabatic. 
However,  heat  transfer  does  take  place  between  the  flows  within  the  mixing 
tube 

7.  The  gas  turbine  provided  constant  mass  flow  of  exhaust  gases  ( m  a=l  89  Ibm/s, 
Ta=910  °F)  at  full  power. 
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B,  THE  IDEAL  ONE-DIMENSIONAL  MODEL 

The  equations  used  for  the  model  come  directly  from  basic  fluid  mechanics 


namely,  Continuity,  Momentum,  Energy  and  State.  Four  control  volumes  were  used  to 

analyze  the  model,  the  first  was  around  the  straight  duct  section,  the  second  around  the 

transition  duct,  the  third  around  the  nozzle  section,  and  the  fourth  around  the  mixing  tube. 

See  Figure  [7]  for  clarification.  It  was  assumed  the  engine  provided  a 

constant  mass  flow 

of  exhaust  gasses  at  a  full  power  setting,  and  that  there  are  no  leaks 

or  sources  of  air  in 

the  ducting  from  the  entrance  (station  “a”)  through  the  nozzle 

Therefore,  from  continuity 

exit  (station  “3”). 

II 

<N 

•s 

II 

(1) 

The  following  equations  were  used  for  the  first  control  volume: 

Continuity 

=  ihj 

(2) 

=  PxAVx 

(3) 

Energy  equation 

P  P 

"  +  "  +Z  =  '  +  '  +Z,+H„ 

P.g  2g  2g  '' 

(4) 

Equations  of  state 

Pa=PaRT, 

(5) 

P,=p,RT, 

(6) 

Bernoulli's  equation 

1  9 

Pta=Pa+l^PaVa 

(V) 

(8) 

Head  loss  due  to  friction 

7-7 

Hn=f  '  “  “ 

'  A.  2g 

(9) 
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Where  the  Darey  frietion  faetor  ( / )  was  assumed  to  be  0.01  for  new  stainless  steel  based 


on  the  duct  hydraulic  diameter,  and  Z.  (i=a,  1 . . .)  is  the  height  of  the  station  cross-section. 


In  similar  fashion,  the  following  equations  were  used  for  the  second  control 
volume: 

Continuity 


m,  =  rh^  (10) 

^2  ~  (11) 


Energy  equation 


P  P 

^  +  -^  +  Z  =^  +  -^  +  Z,+//„ 

p,g  Ig  p,g  2g 


Equation  of  state 


(12) 


Pi  =  PiPTi 


(13) 


Transition  duct  losses 


H 


/2  “ 


(14) 


Where  is  the  total  pressure  loss  coefficient  for  the  transition  duct.  A  detailed 


explanation  of  how  was  determined  is  contained  in  chapter  IV. 


Eigure  6  Plan  view  of  mixing  tube.  Area  A4  is  shaded. 
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m  g ,  Vg  (uniform) 


Figure  7  Station  labeling  of  eduetor  system. 

The  nozzle  seetion  again  used  the  same  equations  as  the  transition  duet  with  the 
exeeption  of  the  total  pressure  loss  eoeffieient.  The  losses  for  the  nozzle  seetion  were 
assumed  to  be  =0.1  sinee  the  aetual  total  pressure  loss  is  a  funetion  of  the  lobed 
nozzle  geometry  as  determined  in  referenee  [5]. 
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The  amount  of  secondary  air  required  was  determined  using  the  following 
equations: 

Equation  of  state 

P,=P,RT,  (15) 

m^=  (16) 

Bernoulli's  equation  for  secondary  air 

P,.=P.+\p,Vi  (17) 

It  was  assumed  that  static  pressures  of  the  primary  flow  at  the  nozzle  exit  and  the 
secondary  flow  at  the  mixing  tube  entrance  were  equal. 

P,=P,  (18) 

Finally,  the  equations  used  for  the  fourth  control  volume  are  as  follows: 

Continuity 

m^+m^=m^  (19) 

=  PeAK  (20) 

4  =  4  +  4  (21) 

For  clarification,  is  shown  as  the  shaded  area  in  Figure  [6].  The  Momentum  equation 

is  given  by 

-  m,V,  -  m,V,  =P,A,+  P,A,  -  P^  (22) 

Energy  equation 

+^r,^)  +  m,(-^RT,  +^Vf)  =  m,(^RT,  +)r;)  (23) 

/-I  2  /-I  2  /-I  2 

Equation  of  state 

Pe  =  PeRTe  (24) 

All  the  equations  previously  mentioned  in  this  section  were  written  into  an 
Engineering  Equation  Solver  (EES)  code  and  solved  simultaneously.  The  EES  software 
solves  the  equations  by  iteration  to  a  specified  tolerance  using  the  known  values  as 
defined  by  the  user,  such  as  ambient  conditions,  gas  turbine  exhaust  mass  flow  rate  and 
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temperature,  and  user  defined  guesses  for  the  unknown  variables.  The  EES  code  that  was 
written  for  this  study  is  contained  in  Appendix  A. 

C.  MASS  FLOW  CORRECTIONS  AND  VALIDATION 

As  was  stated  in  the  assumptions  for  this  study,  the  wall  friction  in  the  mixing 
tube  was  neglected  for  the  ideal  calculations.  This  assumption  results  in  an  overestimate 
of  the  secondary  mass  flow  rate  and  an  underestimate  of  the  mixing  tube  exit 
temperature.  In  order  to  correct  the  model  and  obtain  realistic  secondary  mass  flow  rate 
and  mixing  tube  exit  temperature,  a  correction  factor  (c/)  was  applied  to  the  ideal 
secondary  air  mass  flow  rate  in  the  following  manner 

Korr  =  icf)W  (25) 

Where  the  pumping  ratio  W  is  defined  as 

fh 

_  secondary  (26) 

rii 

primary 

The  equations  for  the  fourth  control  volume  were  then  recalculated  using  the  corrected 
pumping  ratio  to  determine  the  real  mixing  tube  exit  temperature. 

The  correction  factor  was  obtained  from  a  combination  of  empirical  data  found  in 
Keenan  [8]  and  Otis  [6].  The  data  from  Keenan,  shown  in  figure  [8],  are  the  ratio  of 
measured  mass  flow  to  calculated  mass  flow  versus  the  ratio  of  nozzle  total  pressure  to 
ambient  pressure. 


figure  8  Ratio  of  mass  flow  ratios  versus  pressure  ratio,  from  Reference  [6]. 
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A  best- fit  curve  was  applied  to  the  data  to  determine  a  correction  factor.  For  purposes  of 
model  validation,  full  power  for  the  LM2500  was  used,  where  m=153  Ibm/s,  plus  20 
Ibm/s  module  cooling  air,  and  Ta=807  °F,  taken  from  Reference  [7].  Since  the  pressure 
ratio  for  the  DDG  51  eductor  nozzle  is  only  slightly  greater  than  unity,  the  correction 
factor,  from  the  graph  is  0.81.  Based  on  his  model  validation,  Otis  [6]  determined  that 
the  correction  factor  should  be  0.79  in  order  to  produce  the  temperature  found  at  the  exit 
plane  of  the  DDG  5 1  mixing  tube.  Since  the  ideal  model  for  this  study  indicated  that  the 
pressure  ratio  was  also  only  slightly  higher  than  unity,  it  was  decided  to  use  a  correction 
factor  of  0.8.  To  validate  the  model  used  in  this  study,  the  data  used  by  Otis  [6]  for  the 
DDG  51  along  with  the  appropriate  nozzle  loss  coefficient  was  used.  The  exit 
temperature  calculated  by  the  model  differed  by  less  than  1%  compared  to  the  value 
found  in  Reference  [9]  for  DDG  5 1 . 

The  backpressure  that  the  eductor  system  would  develop  was  calculated,  since 
increased  backpressure  decreases  engine  performance.  Engine  backpressure,  as  defined 
in  Reference  [9]  is  the  difference  between  the  static  pressure  measured  at  the  engine 
module  exit  and  ambient  pressure.  Since  the  exhaust  configuration  for  LHA  (R)  from  the 
02  level  to  the  engine  module  is  not  known,  backpressure  is  defined  as  the  difference 
between  the  predicted  static  pressure  at  the  02  level  and  the  ambient  pressure.  The 
backpressure  limit  for  the  LM2500  with  an  eductor  system,  from  Reference  [7]  is  12 
inches  H2O.  For  the  purpose  of  this  study,  the  same  backpressure  limit  was  assumed  for 
the  LM2500+.  Again,  using  information  taken  from  DDG  51,  the  backpressure 
calculated  by  the  model  was  within  18%  of  the  measured  value  from  Reference  [9].  The 
backpressure  calculated  by  Otis  [6]  for  the  DDG  51  was  within  17%  of  the  value  listed  in 
Reference  [9].  Based  on  the  small  error  between  the  calculated  value  and  the  actual  exit 
temperature,  and  the  agreement  in  calculated  backpressures,  it  was  determined  that  the 
code  was  accurate  in  predicting  mixing  tube  exit  temperature  and  secondary  air  mass 
flow. 
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III.  PRELIMINARY  DESIGN  STUDY 


A,  TRADE  STUDY  USING  THE  1-D  CORRECTED  MODEL 

The  corrected  and  validated  eductor  model  discussed  in  the  previous  chapter  was 
used  as  the  primary  design  tool  for  this  study.  The  first  step  in  performing  the  study  was 
to  change  the  values  used  in  the  model  to  those  of  the  LM2500+.  The  values  used  for 
exhaust  mass  flow  {m=  189  Ibm/s  plus  20  Ibm/s  module  cooling  air)  and  turbine  exit 
temperature  (Ta=  910  °F)  were  determined  by  the  Naval  Surface  Warfare  Center, 
Carderock  Division,  these  values  differ  slightly  from  those  quoted  by  General  Electric 
Aircraft  Engines  Division  [16],  where  m  =  189  Ibm/s  and  Ta  =  965  “E. 

The  next  step  was  to  determine  the  available  volume  of  the  existing  EHD  8  uptake 
space.  This  established  the  maximum  size  that  would  be  allowed  for  the  mixing  tube 
since  the  mixing  tube  has  the  largest  cross-section  of  the  exhaust  ducting.  Since  the 
design  for  the  EHA  (R)  superstructure  and  uptake  space  was  not  available  for  this  study, 
it  was  assumed  that  the  uptake  space  for  the  EHA  (R)  was  identical  to  that  of  the  EHD  8. 
With  the  knowledge  of  what  the  maximum  cross-sectional  area  could  be  for  the  mixing 
tube,  several  mixing  tube  cross-sectional  areas  were  chosen  and  the  nozzle  exit  area 
varied  for  each  case.  At  this  point  in  the  design  process,  the  transition  duct  exit  area  was 
neglected  in  order  to  establish  the  nozzle  exit  area.  The  analysis  shown  in  Eigure  [9] 
determined  that  the  mixing  tube  area  to  nozzle  area  (Ag/As)  was  3.1  for  each  case,  where 
the  goal  was  to  have  a  pumping  ratio  of  unity.  It  is  also  evident  in  Eigure  [10]  that  as  the 
mixing  tube  area  decreases  the  backpressure  increases.  Noting  the  trend  of  backpressure 
as  a  function  of  mixing  tube  area,  it  was  decided  to  use  the  largest  possible  mixing  tube 
that  would  fit  into  the  uptake  space. 
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Figure  9  Mixing  tube  exit  area  versus  Nozzle  exit  area  for  a  pumping  ratio  of  unity. 
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Figure  10  Backpressure  versus  mixing  tube  exit  area  for  Ae/A3=3.1  and  W=l. 

Once  the  mixing  tube  area  was  established,  the  next  step  in  the  design  process 
was  to  determine  the  nozzle  exit  area.  Finding  the  nozzle  exit  area  was  performed  in  the 
same  manner,  as  was  the  mixing  tube  area.  By  varying  the  nozzle  exit  area  and  noting 
where  that  area  achieved  a  pumping  ratio  of  unity,  the  correct  nozzle  area  could  be 
determined  for  the  given  primary  flow  conditions.  Figure  [11]  is  a  graphical 
representation  of  both  backpressure  and  exit  temperature  as  a  function  of  nozzle  exit  area. 
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Figure  11  Backpressure  and  exit  temperature  versus  nozzle  exit  area  for  a  constant 
mixing  tube  exit  area. 

Now  that  the  nozzle  exit  area  was  known,  the  exit  area  of  the  transition  duct  could 
be  fixed  and  the  final  design  complete.  The  transition  duct  exit  area  was  set  to  be  larger 
than  the  nozzle  exit  area  so  that  the  cross-section  constantly  decreased  through  the 
nozzle. 

The  selected  nozzle  area  and  mixing  tube  area  resulted  in  an  exit  temperature  that  was 
12%  lower  than  that  achieved  by  the  DDG  51  eductor  system,  at  the  cost  of  a  6  in  H2O 
increase  in  backpressure,  as  compared  to  a  straight  exhaust  duct  without  cooling. 

B.  EDUCTOR  DESIGN 

Now  that  the  general  cross-sectional  areas  have  been  selected,  it  was  necessary  to 
make  them  fit  within  the  available  space  in  the  existing  uptake  space.  This  was 
accomplished  by  squeezing  the  two  diesel  generator  exhausts  into  one-fifth  of  the 
available  area  of  the  uptake  space,  and  flattening  out  the  gas  turbine  exhaust  with  an  oval- 
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to -rectangular  transition  duct.  In  order  to  accomplish  the  rerouting  of  the  diesel  exhaust, 
a  similar  study  was  conducted  using  the  same  1-D  model,  substituting  the  diesel  exhaust 
parameters.  Results  of  the  diesel  exhaust  study  are  contained  in  Appendix  B.  The 
resultant  design  proposed  for  the  LHA  (R)  uptakes  are  shown  in  Figures  [12]  and  [13]. 


GT  EXHAUST 


Figure  12  Plan  view  of  proposed  LHA  (R)  exhaust  system. 
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The  proposed  design  begins  at  the  02  level.  The  gas  turbine  exhaust  eonsists  of 
the  existing  flat  oval  duet,  an  oval-to-reetangular  transition  duet,  a  lobed  nozzle,  and 
mixing  tube.  The  diesel  exhaust  consists  of  the  existing  flat  oval  duct  that  transitions 
through  an  “S”  duct  into  a  flat  rectangle  with  a  lobed  nozzle  and  mixing  tube.  The 
secondary  air  intakes  for  the  gas  turbine  and  diesel  exhausts  are  to  be  separated  so  that 
the  diesel  eductor  is  not  starved  by  the  gas  turbine  eductor. 

C.  BACK  PRESSURE  COMPARISON 

The  following  is  a  comparison  of  backpressure  between  the  proposed  LHA  (R) 
eductor  system  and  an  uncooled  LHA  (R)  gas  turbine  exhaust.  Figure  [14]  shows  the 
change  in  static  and  total  pressure  for  each  section  of  ducting,  measured  in  “H2O.  As 
previously  mentioned  in  Chapter  II,  backpressure  is  defined  as  the  difference  between  the 
predicted  static  pressure  at  the  02  level  and  the  ambient  pressure.  The  addition  of  the 
eductor  components  adds  approximately  6  “H2O  backpressure  to  the  system  as  proposed. 
The  backpressure  in  the  eductor  system  could  be  lowered  at  the  cost  of  increased  exit 
temperature,  see  Figure  [11].  Conversely,  the  exit  temperature  could  be  reduced  further 
if  increased  backpressure  is  allowed. 

D,  THE  LOBED  NOZZLE 

A  lobed  nozzle  was  selected  for  the  eductor  system  because  of  its  enhanced 
mixing  characteristics.  The  lobed  nozzle  takes  advantage  of  differences  in  velocity  and 
density  between  the  primary  and  secondary  flows.  As  the  higher  velocity,  lower  density 
primary  flow  meets  the  slower,  higher  density  secondary  flow,  a  streamwise  vortex  is 
generated,  as  shown  Figure  [5].  The  vortices  result  in  more  rapid  mixing  of  the  two 
flows.  This  vortex  driven  mixing  allows  for  a  reduction  in  the  length  of  the  mixing  tube. 
Reduction  of  the  larger  mixing  tube  will  result  in  space  and  weight  savings  over 
traditional  eductor  systems. 
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Figure  14  Comparison  of  backpressures  for  cooled  and  uncooled  LHA  (R)  exhaust. 
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IV.  TRANSITION  DUCT  DESIGN 


A,  CREATING  THE  TRANSITION 


The  transition  required  going  from  a  flat  oval  eross-seetion  to  a  flat  reetangular 
eross-seetion  with  an  aspeet  ratio  (L/W)  of  approximately  6:1,  with  the  “x”  axis  shifting 
to  one  side,  and  an  inlet  to  outlet  area  ratio  (A1/A2)  of  unity.  Though  many  referenees 
exist  for  eireular-to-reetangular  transitions  (e.g.  Referenee  [10]),  they  do  not  address 
transitions  of  the  eomplexity  required  for  the  duet  used  in  this  study.  Complex  transition 
geometries  are  employed  on  stealth  aireraft.  However,  the  data  for  those  geometries  are 
highly  elassified.  Therefore,  it  was  determined  to  use  the  AR630  transition  duet.  Figure 
[15]  found  in  Patriek  and  MeCormiek  [11],  as  a  starting  point.  The  AR630  duet  eross- 
seetional  areas  are  determined  mathematieally  by  use  of  a  super  ellipse.  The  equation  for 
a  super  ellipse  is 


f-1 

UJ 


+ 


(27) 


Figure  15  AR630  transition  duet,  from  Referenee  [11]. 
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The  challenge  would  be  to  apply  a  super  ellipse  equation  to  begin  as  a  flat  oval  and 
transition  to  a  flat  rectangle  with  a  longitudinal  axis  that  shifts  approximately  2.5  feet. 
The  modified  super  ellipse  equation  used  in  this  study  is 


y  +  c 


\my 


+ 


y  a  J 


=  1 


(28) 


where  the  coefficients  my,  nz,  a,  b  and  c  are  all  functions  of  axial  distance.  The 
coefficients  versus  axial  distance,  along  with  the  code  that  was  used  to  calculate  the 
transition  geometry  are  listed  in  Appendix  C.  Figure  [16]  illustrates  the  geometry 
produced  by  equation  (28)  for  this  study. 


Figure  16  Transition  duct  for  proposed  LHA  (R)  exhaust  system. 

B,  TOTAL  PRESSURE  LOSS  IN  THE  TRANSITION 

To  produce  a  more  accurate  estimate  of  the  performance  for  the  eductor  design, 
the  total  pressure  loss  through  the  transition  duct  must  be  known.  Since  the  losses  for  the 
duct  designed  in  the  previous  section  are  not  known,  it  was  decided  to  use  the  losses 
found  in  the  AR630  duct.  The  geometries  are  similar,  in  that  the  aspect  ratios  and  the 
area  ratios  for  both  ducts  are  nearly  6:1  and  1:1  respectively.  Although  Patrick  and 
McCormick  [11]  provided  a  detailed  report  on  the  AR630  duct,  they  did  not  produce  a 
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total  pressure  loss  eoeffieient  for  the  entire  duet.  Rather,  the  report  showed 
dimensionless  veloeity  and  pressure  profiles  for  both  inlet  and  outlet  eross-seetions.  The 
report  also  ineluded  tabular  data  that  was  used  to  generate  the  outlet  profdes.  The  total 
pressure  loss  eoeffieient  in  the  AR630  duet  was  determined  by  the  equation 


P  -P 

^  t\  ^  tref 


a 


ref 


(29) 


and 


Qref  ~ 


(30) 


where  Vref=100  ft/s  and  Ptref=14.7  psia. 

The  total  pressure  at  the  inlet  of  the  AR630  duet  was  determined  by  generating 
similar  dimensionless  veloeity  and  total  pressure  profiles  in  MATLAB.  In  order  to  obtain 
the  slug  profdes  found  in  Figures  [17]  and  [18],  a  super  ellipse  of  the  form  found  in 
equation  (27),  (where  a=l,  b=l,  m=10  and  n=10)  was  used.  Figures  [18]  and  [19]  are  the 
inlet  veloeity  and  total  pressure  profiles,  for  the  AR630  duet,  generated  by  the  eode. 
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AR630  inlet  velocity  profile 


Figure  19  AR630  inlet  velocity  profile  generated  in  MATLAB. 
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AR630  inlet  total  pressure  profile 
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Figure  20  AR630  inlet  total  pressure  profile  generated  in  MATLAB. 

Onee,  the  profiles  were  caleulated,  the  pressure  was  integrated  over  the  cross- 
section  using  the  equation 

p\\v{x,y)PXx,y)dA 

p^^V{x,y)dA 

The  double  integral  was  calculated  using  MATLAB.  The  results  of  the  calculation 
yielded  the  total  pressure  at  the  duct  inlet. 

Similarly,  the  tabular  data  provided  for  the  outlet  to  the  AR630  duct  was 
integrated  over  the  exit  cross-section  yielding  the  exit  total  pressure.  The  inlet  and  outlet 
total  pressures  were  then  substituted  into  Equation  (29)  to  determine  the  total  pressure 
loss  coefficient  )  for  the  duct,  which  was  estimated  as,  =0.2  .  The  code  used  to 
calculate  the  loss  coefficient  is  provided  in  Appendix  D. 
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C.  VALIDATION  OF  THE  LOSS  COEFFICIENT 

To  ensure  that  integrating  the  inlet  and  outlet  pressure  of  the  duet  in  the  manner 
deseribed  in  the  previous  seetion  produeed  the  eorreet  loss  eoeffieient,  a  simple  model 
was  developed  to  test  the  eode.  The  model  was  a  simple,  gradually  eonverging 
eoneentrie  nozzle  with  a  flow  area  ratio  of  2:1.  From  White  [12],  for  sueh  a  nozzle,  the 
total  pressure  at  the  nozzle  exit  is  95%  that  of  the  total  pressure  at  the  inlet.  Using  the 
super  ellipse  eode  to  develop  the  dimensionless  veloeity  and  pressure  profiles  for  the  inlet 
and  outlet  of  the  nozzle  and  integrating  over  eaeh  eross-seetion,  then  substituting  the  total 
pressures  into  equation  (29)  resulted  in  a  loss  eoeffieient  of  0.049.  The  ealeulated  and 
known  values  for  the  loss  eoeffieient,  for  the  model  ease,  differ  by  3%.  Based  on  the  low 
error,  the  eode  used  to  determine  the  loss  eoeffieient  for  the  AR630  duet  produeed 
aeeurate  results. 
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V.  PLUME  RADIATION  STUDY 


A,  OVERVIEW  AND  ASSUMPTIONS 

The  purpose  of  this  study  is  to  determine  how  much  the  radiant  intensity  of  a  gas 
turbine  exhaust  plume  can  be  reduced  if  the  mixing  tube  exit  geometry  is  a  high  aspect 
ratio  slot  vice  the  current  circular  cross-section.  The  study  will  focus  on  CO2  and  H2O  in 
the  3-5  pm  band,  since  both  CO2  and  H2O  emit  strongly  in  this  band  for  the  temperature 
range  of  the  exhaust  gases.  Several  assumptions  were  made  in  order  to  estimate  the 
plume’s  radiant  intensity,  including  that  the  plume  is  homogeneous  and  isothermal,  and 
does  not  contain  soot.  Additionally,  it  was  assumed  that  the  atmosphere  was  a  non¬ 
participating  medium. 

The  study  used  empirical  data  taken  from  Reference  [9]  for  DDG  51  and  the 
LM2500  gas  turbine  engine  at  full  power  from  Reference  [7].  The  data  from  Reference 
[9]  provided  a  real  plume  temperature,  while  Reference  [7]  provided  the  information 
necessary  to  determine  the  partial  pressures  of  CO2  and  H2O  in  the  DDG  51  exhaust 
plume.  Since  the  LM2500  uses  either  JP-5  or  Diesel  Fuel  Marine,  both  of  which  consist 
of  several  different  types  of  fuels  and  have  no  specific  chemical  designation,  it  was 
assumed  that  the  chemical  composition  of  the  fuel  was  CioHig.igs.  The  composition 
resulted  in  partial  pressures  of  5%  for  both  C02  and  H20  given  a  fuel-to-air  ratio  if)  of 
0.0224. 

The  following  assumptions  were  made  in  order  to  estimate  the  intensity  of  a  gas 
turbine  exhaust  plume; 

1 .  The  plume  was  isothermal  and  homogeneous. 

2.  The  plume  consisted  of  the  by-products  of  complete  combustion  only. 

3 .  The  atmosphere  was  a  non-participating  medium. 

4.  The  chemical  composition  of  the  fuel  used  in  U.  S.  Navy  marine  gas  turbines 
was  C10H19.185. 
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B, 


THE  EXPONENTIAL  WIDE  BAND  MODEL 


According  to  Modest  [13],  the  Exponential  Wide  Band  Model  “is  by  far  the  most 
sueeessful  of  the  wide  band  models”,  with  errors  of  up  to  +/-  20%  when  eompared  with 
empirieal  data.  Therefore,  the  Exponential  Wide  Band  Model  will  be  used  to  ealeulate 
the  emissivity  of  the  plume  for  use  in  the  equation: 

E  =  £aT^  [wW-sr]  (32) 

where  E  is  the  emissive  power  of  the  plume,  s  is  the  emissivity,  a  is  the  Stefan- 
Boltzmann  eonstant  and  Te  is  the  plume  temperature.  The  radiant  intensity  of  the  plume 
ean  be  determined  from 

/  =  —  [wW-sr]  (33) 

n 

where  1  is  the  radiant  intensity  of  the  plume.  The  equations  used  to  ealeulate  emissivity 
are  as  follows: 

The  band  width  parameter 


The  band  overlap  parameter 


(35) 


where  the  term  in  parenthesis  is  taken  from  Eigures  [21]  and  [22],  for  H2O  and  CO2 
respeetively. 
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Figure  21  Temperature  Dependenee  of  the  line  overlap  parameter,  y,  and  band  strength 
parameter,  a,  for  water  vapor,  from  Reference  [13]. 
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Normalized  band  strength  parameter  uIOq 


Figure  22  Temperature  dependenee  of  the  line  overlap  parameter,  y,  and  band  strength 
parameter,  a,  for  earbon  dioxide,  from  Referenee  [13]. 

The  effeetive  pressure  is  given  by 


P,  =  —  l  +  (h-l)^ 


where  po  is  1  atm,  pa  is  the  total  pressure  of  the  plume  and  p  is  the  partial  pressure  of  the 
gas,  either  CO2  or  H2O.  The  line  overlap  parameter,  P  is  given  by 

P  =  YP.  (37) 


The  pressure  path  length  is 


X  =  P.s 


where  s  is  the  path  length  through  the  plume  as  shown  in  Figure  [24]  and  pa  is  the  density 
of  the  absorbing  gas.  The  gas  band  optieal  thiekness  is  given  by 


T  -  ° 

^  o 
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The  non-dimensional  band  absorbtanee,  A*,  ean  now  be  determined  from  the  following 
table: 

TABLE  9.2 

Exponential  wide  band  correlation  for  an  isothermal  gas _ 


)3  <  1 

0  <  To  <  i3 

II 

« 

Linear  regime 

/3  <  -Jb  <  l/)3 

A'  =  2x/7b/3  -  /3 

Square  root  regime 

l/j3  <  7t)  <  00 

A*  =  ln(ib]8)  +  2  - 

Logarithmic  regime 

1 

0  <  1b  ^  1 

A*  =  7b 

Linear  regime 

1  ^  H)  <  00 

A*  =  In  7b  +  1 

Logarithmic  regime 

a,  and 

w  from  Table  9.3  and  equations  (9.69)  through  (9.76),  Tq  =  aXito. 

Figure  23  Exponential  wide  band  eorrelation  for  an  isothermal  gas,  from  Referenee  [13]. 


For  eaeh  of  the  preeeding  equations,  the  eorrelation  parameters  of  n,  b,  Uo,  Yo  and  cOo  are 
all  taken  from  table  9.3  of  referenee  [13]. 

The  total  band  absorbtanee  is  then  ealeulated  from 

A  =  A*0)  (40) 

Finally,  the  emissivity  is  given  by 


^ A^ 


"bt] 


Where 


^br, 


is  found  by  interpolation  of 


(41) 


in  Appendix  C  of  referenee  [12]. 


The  preeeding  equations  were  written  into  a  MATFAB  eode  to  assist  in  making 
the  ealeulations.  The  eode  was  verified  by  using  Example  9-6  in  Modest  [13].  The 
values  ealeulated  in  by  the  eode  were  identieal  to  those  found  in  Example  9-6  therefore, 
the  eode  was  used  for  the  study.  A  eopy  of  the  eode  is  eontained  in  Appendix  E. 


39 


C.  PLUME  INTESITY  VERSUS  PATH  LENGTH 


Using  the  code  described  in  the  previous  section,  the  plume’s  radiant  intensity 
was  estimated  for  various  aspect  ratios  of  the  same  cross-sectional  area.  The  first  step 
was  to  determine  the  fuel-to-air  ratio,  for  this  study  performance  data  for  the  LM2500  at 
full  power  was  used.  Based  on  an  exhaust  mass  flow  rate  of  rh  exhaust  =153  Ibm/s  and 
specific  fuel  consumption,  SFC,  of  0.37  Ibm/hp-hr,  the  fuel-to-air  ratio  could  be 
determined  from  the  following  equations: 

/  =  (42) 


^fuei  -  *  Ps  where  Ps  =  32,000  hp 

^fuel  ^air  ~  ^exhaust 


(43) 

(44) 


The  resultant  fuel-to-air  ratio  was  0.0224. 

The  next  step  was  to  determine  the  proper  stoichiometry  for  complete 
combustion.  Assuming  there  was  no  excess  air  used  in  the  combustion  process,  the 
balanced  reaction  would  be 


C10//19 185  +I4.8O2  +  67.7^2-  >  IOCO2  +9.59/72(9  +  2402  +  67.7^2  (45) 

Since  the  LM2500  requires  44.6  Ibm  of  air  for  every  Ibm  of  fuel  burned,  and  assuming  air 
is  20%  O2  and  80%  N2,  the  amount  of  excess  air  used  in  the  actual  combustion  process 
can  be  determined.  The  excess  air  can  then  be  added  to  the  combustion  process  in 
Equation  (45)  resulting  in  the  following  stoichiometric  reaction 

C10//19 185  +38.8O2  +177.5^2-  >  IOCO2  +  9.59//2O  +  24O2  +177.5^2  (46) 

Where  the  partial  pressures  of  the  exhaust  gases  can  be  determined  from  the  mole 
fractions  as  follows 


#  moles  „ 


#  moles 


-{\atm) 


(46) 


total 


The  partial  pressures  of  CO2  and  H2O  are  5%  and  5%  respectively.  The  partial  pressures 
and  temperature  of  the  plume  were  then  entered  into  the  code  to  determine  the  radiant 
intensity  of  the  plume. 
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The  first  path  length,  s,  was  that  of  the  eurrent  DDG  5 1  geometry  where  the  path 
length  was  through  the  plume  axis.  The  eireular  eross-seetion  was  then  eonverted  to  a 
square  of  the  same  area.  Next  the  aspeet  ratio  (AR)  was  varied  from  1:1  to  40:1.  The 
minor  axis  was  the  path  length  of  interest,  Figure  [25],  as  it  was  assumed  that  the  exhaust 
system  would  be  oriented  in  sueh  a  way  as  to  present  the  minimum  path  length  to  an 
ineoming  missile. 


Figure  25  Plan  view  illustrating  exhaust  orientation  to  IR  Seeker. 
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The  study  revealed,  as  expected,  that  the  radiant  intensity  of  the  plume  decreases 
as  path  length  decreases.  Figure  [26]  illustrates  the  dependence  of  plume  intensity  on 
path  length  for  a  cooled  and  uncooled  plume,  while  Figures  [27]  and  [28]  show  the 
emissivities  of  the  individual  gas  at  the  specified  wavelength.  Based  on  the  following 
graphs,  the  radiant  intensity  of  the  cooled  exhaust  can  be  reduced  to  one-eighth  that  of  the 
uncooled  plume  for  the  circular  cross-section.  If  the  total  gas  emissivities  for  the  cooled 
and  uncooled  plume  were  equal,  it  would  be  expected  that  I  cooled  =  l/16/„„coo/erf  since  the 
temperature  of  the  cooled  exhaust  is  nearly  half  that  of  the  uncooled,  as  given  by 
Tuncoded  =  ^T'^ooied  ^nd  I ■  Howcvcr,  the  total  gas  emissivities  are  not 
equal.  Figures  [27  d]  and  [28  d]  show  that  the  total  gas  emissivity  of  the  uncooled  plume 
is  slightly  more  than  double  that  of  the  cooled  plume,  so  that  Icooied  =  ^/^luncooied- 


Figure  26  Total  radiation  intensity  versus  path  length  of  a  plume. 
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Figure  27  Gas  emissivity  versus  path  length  for  each  spectral  band  (uncooled  plume  at 
5%  CO2  and  5%  H2O). 
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Figure  28  Gas  emissivity  versus  path  length  for  each  spectral  band  (cooled  plume  at  5% 
CO2  and  5%  FI2O). 

The  preceding  graphs  show  that  cooling  the  plume  significantly  reduces  its 
infrared  signature,  and  converting  the  exit  geometry  from  a  circular  cross-section  to  a 
long  thin  rectangle  can  reduce  the  radiant  intensity  of  the  plume.  Figure  [29]  shows  that 
for  a  cooled  plume,  the  radiant  intensity  can  be  reduced  by  40%  if  a  40:1  aspect  ratio  slot 
were  used  in  place  of  a  circular  cross-section,  this  would  result  in  a  slot  that  is 
approximately  40  ft  by  1  ft  for  DDG  51.  In  comparison,  a  more  reasonable  aspect  ratio 
of  6:1  results  in  a  20%  reduction  in  plume  intensity.  Again,  the  results  of  this  study  are 
based  on  the  assumption  that  the  plume  is  isothermal  and  homogeneous.  Based  on  the 
results  of  the  study,  the  shipbuilder  must  be  willing  to  compromise  between  the  possible 
loss  of  topside  space  and  weight  gain  in  order  to  reduce  the  radiant  intensity  of  the  plume 
a  significant  amount. 


H20  at  2.7  micron  C02  at  4.3  micron 
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Figure  29  Total  gas  intensity  versus  path  length  for  cooled  plume  showing  various  aspect 
ratios  and  corresponding  intensities. 
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VI.  SUMMARY  AND  CONCLUSIONS 


A,  SUMMARY 

A  one-dimensional  model  of  an  enhaneed  mixing  eduetor  system  was  developed 
and  used  to  prediet  exhaust  temperature,  seeondary  air  mass  flow  rate  and  baekpressure. 
The  predietions  obtained  from  the  model  were  then  used  to  determine  the  preliminary 
design  of  a  gas  turbine  exhaust  system  for  the  LHA  (R)  program.  The  design  ineluded  an 
oval-to-reetangular  transition  duet.  The  geometry  of  the  transition  was  determined  using 
a  super-ellipse  and  the  total  pressure  losses  within  the  transition  were  determined  from 
empirieal  data  taken  from  a  similar  duet  design.  The  duet  losses  were  ineluded  in  the 
exhaust  system  model  to  enhanee  the  predietions.  A  lobed  nozzle  was  seleeted  for  the 
design  in  order  to  reduee  the  mixing  tube  length.  Additionally,  a  eode  was  written,  that 
used  an  Exponential  Wide  Band  Model  to  study  the  effeet  of  path  length  on  exhaust 
plume  radiant  intensity,  and  showed  that  ehanging  the  aspeet  ratio  of  the  mixing  tube  exit 
plane  ean  further  reduee  the  plume’s  infrared  signature. 

B,  CONCLUSIONS 

The  following  eonelusions  were  obtained; 

1 .  An  enhaneed  mixing  eduetor  system  that  fits  within  the  existing  framework  of 
the  LHA  (R)  uptake  spaee  was  designed.  The  proposed  system  provides  a 
46%  reduetion  in  exhaust  temperature,  an  improvement  of  10%  eompared  to 
DDG51. 

2.  An  estimate  of  the  proposed  system  baekpressure  was  ealeulated  and 
eompared  with  an  uneooled  system.  The  inerease  in  baekpressure  developed 
by  the  eduetor  system  was  approximately  6  “H2O  whieh  is  within  the  limit  of 
the  engine.  Additionally,  an  oval-to-reetangular  transition  geometry  was 
provided.  The  estimated  total  pressure  loss  eoeffieient  for  the  transition  was 
estimated  to  be  0.2. 
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3.  The  radiant  intensity  of  the  plume  can  be  reduced  greater  than  40%  for  slotted 
mixing  tube  geometries  with  an  aspect  ratio  of  greater  than  40:1.  However, 
due  to  space  and  weight  constraints  a  more  realistic  reduction  in  signature  of 
up  to  20%  can  be  achieved  for  an  aspect  ratio  of  6: 1 . 
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APPENDIX  A 

1-D  CORRECTED  MODEE 

{Thesis  Project:  EHA(R)  Exhaust  Design  Study} 

{Written  by  Steve  Dudar,  MSME  Dec  03} 

{See  Figure  [7]  for  clarification  of  symbology} 

{CONSTANTS} 

{Gas  constant  R  [J/kg-K],  assumed  exhaust  gases  were  air} 
R=287 

{Specific  heat  ratio  (assumed  constant),  gamma=Cp/Cv} 
gamma=1.4 

{Acceleration  due  to  gravity  [m/s^2]} 
g=9.81 

{Density  of  fresh  water  [kg/m^3]} 
rho_h2o=1000 

{Assumed  ISO-STP  at  sea  level,  T  [K],  P  [Pa]} 

T_amb=288 

P_amb=101325 

{P_4=P_amb} 

{Friction  factor,  for  new  stainless  steel} 
f=0.01 

{PRIMARY  FEOW  at  entrance} 

{From  EM-2500+  technical  data  at  given  conditions} 

{CONTROE  VOEUME  1} 

{Exhaust  mass  flow  at  turbine  exit  [kg/s]} 

m_dot_a=94.8 

m_dot_a=rho_a*A_a*V_a 

{Exhaust  temperature  at  turbine  exit  [K]} 

T_a=760 

{Stack  cross-sectional  area  [m^2]} 

A_a=3.6 

{A_l=A_a} 

{Hydraulic  diameter  [m]} 

D  ha=1.95 
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{Height  of  section  1  [m]} 

Z_a=24.6 

Z_l=36.3 

{Equations  of  state} 

P_a=rho_a*R*T_a 

P_l=rho_a*R*T_l 

{Continuity} 
m_dot_a=m_dot_  1 
m_dot_l=rho_a*A_l  *V_1 

{Bernoulli's  equation  at  entrance} 

P_ta=P_a+0.5  *rho_a*  V_a^2 
P_tl=P_l+0.5*rho_a*V_1^2 

{Energy  equation} 

P_a/ (rho_a*  g)+(  V_a^2)/ (g  *  2)+Z_a=P_  1  /  (rho_a  *  g)+(  V_  1  ^2)/(2  *  g)+Z_  1  +H_f  1 
H_fl=P((Z_l-Z_a)/D_ha)*(V_a^2)/(2*g) 

H_lair=H_fl 
SG_1  =rho_a/rho_h2o 

H_lh2o=H_lair*SG_l/0.3048*12  {[inches  H20]} 

{CONTROE  VOEUME  2} 

{Cross  sectional  area  [m^2]  at  transition  exit  plane} 

A_2=3.37 

{Height  of  section  2} 

Z_2=39.32 

{Equation  of  state} 

P_2=rho_a*R*T_2 

{Continuity} 
m_dot_2=m_dot_a 
m_dot_2=rho_a*  A_2  *  V_2 

{Energy  equation} 

P_l/(rho_a*g)+(V_1^2)/(g*2)+Z_l=P_2/(rho_a*g)+(V_2^2)/(2*g)+Z_2+H_f2 


zeta_t=0.2  {Minor  loss  coefficient  for  transition  duct} 

H_f2=zeta_t*(V_2^2)/(2*g) 

H_2air=H_f2 

H_2h2o=H_2air*SG_l/0.3048*12  {[inches  H20]} 
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{CONTROL  VOLUME  3} 

{Cross  sectional  area  at  nozzle  exit  plane} 

{A_3=2.276  }  {VARIABLE} 

{Height  of  section  3} 

Z_3=40.54 

{Equation  of  state} 

P_3=rho_a*R*T_3 

{Continuity} 
m_dot_3 =m_dot_a 
m_dot_3=rho_a*  A_3  *  V_3 

{Energy  equation} 

P_2/(rho_a*g)+(V_2^2)/(g*2)+Z_2=P_3/(rho_a*g)+(V_3^2)/(2*g)+Z_3+H_f3 

zeta_n=0. 1  {Minor  loss  coefficient  for  nozzle} 

H_B=zeta_n*((V_3^2)/(2*g)) 

H_3air=H_f3+((V_3^2)/(2*g)) 

H_3h2o=H_3air*SG_l/0.3048*12  {[inches  H20]} 

A_l=A_a 

{SECONDARY  EEOW} 

{Assume  that  secondary  flow  is  air  at  ISO-STP} 

T_4=T_amb 

P_t4=P_amb 

P_4=P_3 

{Secondary  air  mass  flow  [kg/s]} 
m_dot_4=rho_4  *  A_4  *  V_4 

{BERNOUEEI'S  EQUATION  AT  THE  MIXING  TUBE  INLET} 
P_t4=P_4+(0.5  *rho_4*(V_4^2)) 

{Equation  of  state} 

P_4=rho_4  *  R*  T_4 

{TERTIARY  EEOW} 

{Assume  Pe=Pamb} 

P_e=P_amb 

m_dot_e=rho_e  *  A_e  *  V_e 
A_e=(A_3+A_4) 

P  e=rho  e*R*T  e 
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A_e=7 

{continuity} 

m_dot_a+m_dot_4=m_dot_e 
{Momentum  equation} 

m_dot_e*  V_e-m_dot_3  *  V_3-m_dot_4*  V_4=P_3  *  A_3+P_4*  A_4-P_e*  A_e 
{Energy  equation} 

m_dot_3  *((gamma/(gamma-  l))*R*T_3+0.5  *V_3^2)+m_dot_4*((gamma/(gamma- 
1 ))  T_4+0 . 5  *  V_4^2)=m_dot_e*  ((gamma/(gamma- 1 ))  *R*  T_e+0 .5  *  V_e^2) 

{STATIC  BACK  PRESSURES} 

DEETAP  =P_ta-P_3 

{Pumping  ratio} 

W=m_dot_4/ m_dot_a 

{Mass  flow  eorreetion} 
m_dot_4_real=0 . 8  *  m_do  t_4 
m_dot_4_real=rho_4  *  A_4  *  V_4_real 
W_real=m_dot_4_real/ m_dot_a 
m_dot_e_real=m_dot_a+m_dot_4_real 
m_dot_e_real=rho_e*A_e*V_e_real 

{Momentum  equation} 

m_dot_e_real*V_e_real-m_dot_3*V_3- 

m_dot_4_real*V_4_real=P_3*A_3+P_4_real*A_4-P_e*A_e 

{Energy  equation} 

m_dot_3  *((gamma/(gamma-  l))*R*T_3+0.5  *  V_a^2)+m_dot_4_real*((gamma/(gamma- 
1  ))*R*T_4+0.5  *  V_4_reaU2)=m_dot_e_real*((gamma/(gamma- 
1  ))*R*T_e_real+0.5  *  V_e_reaU2) 


54 


APPENDIX  B 

DIESEL  EXHAUST  DESIGN  STUDY 

Since  the  diesel  exhaust  had  to  be  rearranged  in  order  to  aoeommodate  the 
proposed  gas  turbine  exhaust  eduetor  system,  it  was  deeided  to  propose  a  similar  eduetor 
system  for  the  diesel  exhaust.  Following  the  same  methodology  as  deseribed  in  ehapter 
III,  and  using  the  parameters  for  a  Colt-Pielstiek  PA6  diesel  engine  operating  at  full 
power,  nozzle  exit  area  was  determined.  It  was  not  neeessary  to  determine  the  mixing 
tube  exit  area  using  the  eode  sinee  the  area  was  eonstrained  by  the  gas  turbine  exhaust 
mixing  tube.  Figure  [25]  shows  baekpressure  and  exit  temperature  versus  nozzle  exit 
area.  In  order  to  determine  the  baekpressure  plaeed  on  the  engine  due  to  the  eduetor 
system,  a  nozzle  loss  eoeffieient  and  “S”  duet  loss  eoeffieient  of  0.1  and  0.2,  respeetively 
was  assumed.  The  resultant  exit  temperature  was  redueed  by  30%  with  an  inerease  of  2  in 
H2O,  as  eompared  to  the  uneooled  diesel  exhaust. 


Figure  30  Baekpressure  and  exit  temperature  versus  nozzle  exit  area  (SSDG). 
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APPENDIX  C 

TRANSITION  DUCT  DESIGN 

Coefficients  for  transition  duct  super-ellipse,  equation  (28) 


X[ft] 

a  [ft] 

b[ft] 

0.0000 

4.7900 

2.2500 

0.5000 

4.9130 

2.1938 

1.0000 

5.0360 

2.1375 

1.5000 

5.1590 

2.0813 

2.0000 

5.2820 

2.0250 

2.5000 

5.4050 

1.9688 

3.0000 

5.5280 

1.9125 

3.5000 

5.6510 

1.8563 

4.0000 

5.7740 

1.8000 

4.5000 

5.8970 

1.7438 

5.0000 

6.0200 

1.6875 

5.5000 

6.1430 

1.6313 

6.0000 

6.2660 

1.5750 

6.5000 

6.3890 

1.5188 

7.0000 

6.5120 

1.4625 

7.5000 

6.6350 

1.4063 

8.0000 

6.7580 

1.3500 

8.5000 

6.8810 

1.2938 

9.0000 

7.0040 

1.2375 

9.5000 

7.1270 

1.1813 

10.0000 

7.2500 

1.1250 

c[ft] 

nz 

my 

0.0000 

3.0000 

0.3333 

0.1230 

3.1250 

0.3200 

0.2460 

3.2500 

0.3077 

0.3690 

3.2500 

0.3077 

0.4920 

3.5000 

0.2857 

0.6150 

3.7500 

0.2667 

0.7380 

4.0000 

0.2500 

0.8610 

4.2500 

0.2353 

0.9840 

4.5000 

0.2222 

1.1070 

5.0000 

0.2000 

1.2300 

5.5000 

0.1818 

1.3530 

6.0000 

0.1667 

1.4760 

7.0000 

0.1429 

1.5990 

8.0000 

0.1250 

1.7220 

9.0000 

0.1111 

1.8450 

10.0000 

0.1000 

1.9680 

10.0000 

0.1000 

2.0910 

10.0000 

0.1000 

2.2140 

10.0000 

0.1000 

2.3370 

10.0000 

0.1000 

2.4600 

10.0000 

0.1000 
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Transition  Duct  MATLAB  code: 

%  This  code  ealeulates  the  geometry  of  an  oval-to-reetangular  transition  duet. 

%  The  eode  only  ealeulates  one  half  the  geometry  about  the  plane  of  symmetry  (x  axis). 
%  Written  by  S.  W.  Dudar 

x=zeros(21); 
for  i=2:22 
x=x+0.5; 

z=0:b(i-l)/20:b(i-l); 

y=a(i- 1  )*((  1  -(z.^nz(i- 1  ))/b(i- 1  )^nz(i- 1  )).^ny(i- 1  ))+o(i- 1 ); 
yl=-y+2*o(i-l); 
plot3(x,y,z,x,y,z,x,y,-z,x,y,-z); 
hold  on 
end 
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APPENDIX  D 

TRANSITION  DUCT  TOTAL  PRESSURE  EOSS  COEFFICIENT  MATEAB  CODE 
%  This  MATEAB  code  was  used  to  determine  the  Total  Pressure  loss  eoeffieient  in  the 


%  AR630  transition  duet,  Ref  [10]. 

%  Written  by  S.  W.  Dudar 
elear  all 

%  Create  AR630  inlet  veloeity  profile 

ra=l; 

rb=l; 


nz=10; 

ny=l/nz; 

dr=0;rb/500:rb; 

r_R=dr; 

Ul_Uref=ra.*((l-(r_R.^nz)./rb.^nz).^ny); 

r_R2=-r_R; 

dTheta=0:2*pi/500;2*pi; 

%  Determine  inlet  Total  Pressure 
Pref=14.7; 

Uref=100; 
rho=0. 06762; 

Qref=0.5*rho/32.2*UreU2; 

Ul=Ul_UrePUref; 

kl=Ul_Uref; 

Pt  1  =k  1 .  *  Qref/ 1 44+Pref; 

Pml=0; 


%{[Psia]} 

%{[ft/s]} 

%{[lbm/fU3]} 

%{[Psia]} 


%{[Psia]} 


mdotl=0; 
for  i=l;501 
for  j=l:501 

Pml=Ptl(j)*Ul(j)*  r_R(i)  *  dr(j)  *  dTheta(i)+Pm  1 ; 
mdotl=U  1  (j)*r_R(0*dr(j)*dTheta(i)+mdotl ; 
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end 

end 

PTl=Pml/mdotl 
%  Determine  exit  Total  Pressure 
k2=xlsread('loss  dataVdP_Qref); 
U2_Uref=xlsread('loss  dataVu_Uref); 
dy=xlsread('loss  dataVdy')'; 
dz=xlsread('loss  data'/dz')'; 

Pt2=k2 .  *  Qref/ 1 44+Pref; 

U2=U2_UrePUref; 

Pm2=0; 
mdot2=0; 
for  i=l;l  1 
for  j=l;ll 

Pm2=Pt2(i,j)*U2(i,j)*dy(i)*dz(j)+Pm2; 

mdot2=U2(i,j)*dy(i)*dz(j)+mdot2; 

end 

end 

PT2=Pm2/mdot2 

%  Determine  Total  Pressure  Loss  Coeffieient 
K=144*(PT1-PT2)/Qref 
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APPENDIX  E 

EXPONENTIAL  WIDE  BAND  MODEL  CODE 


%  {CONSTANTS} 

T_0=100; 

p_0=l; 

p=l; 

sigma=5.67e-8; 


%{[K]} 

%{[atm]} 

%{[atm]} 

%{Stefan-Boltzmann  constant} 


Y=xlsread('Black  Body  Emissivity','Sheet2'); 
x=xlsread('Black  Body  EmissivityVSheetl'); 


%  {VARIABEES} 
T_e=  {user  input} ; 
p={user  input}; 
s=  {variable}; 


%{Exit  Temperature  [K]} 
%{Gas  partial  Pressure  [atm]} 
%{Path  length  [m]} 


%  {Prom  Table  9.3  Modest  [12]} 

eta=; 

n=; 

b=; 

alpha_0=; 

gamma_0=; 

omega_0=; 

E=eta/T_e; 

omega=omega_0  *  sqrt(T_e/T_0); 
gamma=gamma_0  *  2 ; 

P_e=((p/p_0)  *  ( 1  +(b- 1  )*  (p_a/p)))^n; 
beta=gamma*P_e 

rho_a=44*(1000/22.4)*(273/T_e)*p; 

for  i=l:  length} s) 

X(i)=rho_a*s(i); 
tau(i)=alpha_0  *X(i)/ omega 


%{[l/em]} 

%{[m^2/(em*g)]} 

%{[l/em]} 

%{[l/(em*K)]} 


while  beta  <=  1 

if  (tau(i)  <=  beta)  &  (tau(i)  >=  0) 

Al(i)  =  tau(i); 
break 

elseif  (tau(i)  <=  1/beta)  &  (tau(i)  >=  beta) 
Al(i)  =  2*sqrt(tau(i)*beta2)-beta2; 
break 
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elseif  (tau(i)  <  le6)  &  (tau(i)  >=  1/beta) 
Al(i)  =  log(tau(i)*beta)+2-beta2; 
break 
end 
break 
end 

while  beta  >=  1 

if  (tau_co2(i)  <=  1)  &  (tau(i)  >=  0) 

Al(i)  =  tau(i); 
break 

elseif  (tau_eo2(i)  <  le6)  &  (tau(i)  >=  1) 
Al(i)  =  log(tau(i))+l; 
break 
end 
break 
end 

A=Al(i) 

A(i)  =  Al(i)*omega; 

Mi  =  interpl(Y,x,L); 

J2(i)  =  Mi  *A(i); 

epsilon(i)  =  J2(i)/(sigma*T_e); 

1  (i)=(epsilon  (i)*  sigma*  T_e^4)/pi; 


end 

hold  on 
plot(s, epsilon) 
%plot(s,l) 
%plot(s,tau) 
hold  off 
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